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Abstract 

We illustrate how classical chaotic dynamics influences the quantum 
properties at mesoscopic scales. As a model case we study semiclassi- 
cally coherent transport through ballistic mesoscopic systems within the 
Landauer formalism beyond the so-called diagonal approximation, i. e. by 
incorporating classical action correlations. In this context we review and 
explain the two main trajectory-based methods developed for calculat- 
ing quantum corrections: the configuration space approach and the phase 
space approach that can be regarded as less illustrative but more general 
than the first one. 

1 Introduction: quantum transport 
through chaotic conductors 

Among the most important properties characterizing an electronic nanosystem is 
its electrical conductance behavior. Hence gaining knowledge on charge trans- 
port mechanisms, in particular when shrinking conductors from macroscopic 
sizes down to molecular-sized wires or atomic point contacts, has been in the 
focus of experimental and theoretical research throughout the last decade. Such 
a reduction in size and spatial dimensionality goes along with a crossover from 
charge flow in the macroscopic bulk, well described by Ohm's law, to distinct 
quantum effects in the limit of microscopic or atomistic wires. Nanoconductors 
in the crossover regime, often referred to as mesoscopic, frequently exhibit a co- 
existence of both classical remnants of bulk features, combined with signatures 
from wave interference. Such quantum effects usually require low temperatures 
where coherence of the electronic wavefunctions is retained up to micron scales. 
This has lead to the observation of various quantum interference phenomena 
such as quantized steps in the point contact conductance, the Aharonov-Bohm 
effect and universal conductance fluctuations, to name only a few. 

Nonlinear effects can enter into transport through mesoscopic or nano-systems 
twofold: First, as nonlinear I-V characteristics and charge flow far from equilib- 
rium for large enough voltage applied. Second, in the limit of linear response to 
an applied electric field, the intrinsic nonlinear classical dynamics of the unper- 
turbed conductor can govern its transport properties. In this chapter we focus 
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Figure 1: Quantum mechanical wave packet launched into a mesoscopic cavity 
with the geometry of a " dcsymmctrizcd diamond billiard" . The wave packet 
evolution is monitored at times t = 1, 2, 3, 4, and 25, in units of the average time 
between collisions with the walls of a corresponding classical particle (Courtesy 
of A. Goussev). 

on the latter case which is particularly interesting for mesoscopic conductors 
because the nonlinear charge carrier dynamics can influence both the classical 
and, in a more subtle way, the quantum transport phenomena. 

Initially, disordered metals with underlying diffusive charge carrier motion 
were in the focus of interest in mesoscopic matter. Here, we instead address bal- 
listic nano- or mesoscopic conductors where impurity scattering is suppressed. 
The most prominent ballistic systems are nanostructures built from high- mobility 
semiconductor heterostructures, where electrons are confined to two-dimensional, 
billiard-type cavities of controllable geometry; however, such systems are also re- 
alized as atom optics billiards or through wave scattering as optical, microwave, 
or acoustic mesoscopic resonators pQ. 

The quantum dynamics of a single particle in such a geometry is illustrated 
in Fig. Q] showing snapshots of a wave packet after multiples of the average clas- 
sical time between bounces off the billiard walls. The quantum evolution in such 
a mesoscopic geometry, with corresponding chaotic classical dynamics, is char- 
acterized by two main features: The rapid transition from wave packet motion 
following roughly the path of a classical particle to random wave interference at 
larger times and, second, the emergence of wave functions of complex morphol- 
ogy with wave lengths much shorter than the system size. The latter can be 
used to further specify mesoscopic matter, i.e. quantum coherent systems where 
the smallest (quantum) length scale, the de Broglie wave length or Fermi wave 
length Af in electronic conductors, is much smaller than the system size C, i.e. 
l/(k-p£) is a small parameter, in terms of the Fermi momentum fcp = 27t/Af, 
but not fully negligible as in the case of macroscopic systems. 

Such ballistic mesoscopic systems are ideal tools to study the connection 
between (chaotic) classical dynamics and wave interference. Semiclassical tech- 
niques provide this link presumably in the most direct way. Modern semiclassical 
theory is based on trace formulas, sums over Fourier-type components associated 
with classical trajectories. Analogous to the famous Gutzwiller trace formula 
for the density of states [5] , there exist corresponding expressions for quantum 
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transport in the linear response regime. There, semiclassical expressions for the 
conductance have been obtained within the framework of the Landauer-Buttiker 
approach, relating conductance to quantum transmission in nanostructures. 

After early, pioneering semiclassical work by Miller [5J for molecular reac- 
tions and later by Bliimel and Smilansky 4J for quantum chaotic scattering, 
major advances were made in the context of mesoscopic conductance in the 
early nineties by Baranger, Jalabert and Stone [51 All these semiclassical 
approaches were based on, and limited by, the so-called diagonal approxima- 
tion, see below. While most of the features of experimental and numerical 
magneto-conductance profiles could be well explained qualitatively on the level 
of the diagonal approximation, it contained the major drawback that it was 
not current-conserving and hence failed to give correct quantitative predictions 
for the quantum transmission. This was cured about 10 years later when an 
approach was devised to account for off-diagonal contributions to the semiclas- 
sical conductance [7] , thereby achieving unitarity, reflected in (average) current 
conservation, and furthermore agreement with existing predictions from random 
matrix theory (RMT). The applicability of RMT to closed chaotic mesoscopic 
systems was conjectured after numerical simulations in [8]. So it was expected 
to be applicable also to open systems. 

Semiclassical ballistic transport on the level of the diagonal approximation 
was reviewed in detail in Refs. [9j QUI Qj]. In this chapter we hence focus on 
recent progress beyond the diagonal approximation during the last years. This 
serves as a model case which illustrates how chaotic nonlinear dynamics can 
govern quantum properties at nanoscales. 

2 Semiclassical limit of the Landauer transport 
approach 

The Landauer formalism [12j . providing a link between the quantum transmis- 
sion and conductance, has proven to be an appropriate framework to address 
phase-coherent transport through nanosystems. Consider a sample attached to 
two leads of width W% and W2 that support N\ , respectively, N2 current carry- 
ing transverse modes at (Fermi-)energy Ep. For such a two-terminal setup the 
conductance reads at very low temperatures [12j 



with N\ — W\\phnEp j (hw) and an analogous relation for N2. Here, g s = 
2 accounts for spin degeneracy, and the t nm (E) are transmission amplitudes 
between incoming channels m and outgoing channels n in the leads at energy 
E. They can be expressed in terms of the projections of the Green function of 
the scattering region onto the transverse modes <j) n (y') and 4> m (y) in the two 
leads [15] : 



t„ m (E) = -i%„« m ) 1/2 / Ay' / &y<t>* n {y')<t> m {y) G(x',y',x,y;E) . (2) 



Here x and x', respectively, denote the direction along the leads and v„, v rn the 
corresponding longitudinal velocities. The integrals in Eq. ^ are taken over 



G(E F ) = g s e -T(E F ) =5sy EE I W#f)| 2 



2 Ni N 2 



(1) 
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the cross sections of the (straight) leads at the entrance and the exit. 




Ep I eV channel number 

Figure 2: Total quantum transmission (as a function of energy Ep, on the 
right in units of the channel number in the leads of width W\ = W2 — W) 
for transport through a phase-coherent graphene-based quantum dot, see inset. 
The fluctuating line in the left panel is the full quantum transmission at zero 
magnetic field. In the right panel, the straight solid and dashed line denote the 
averaged transmission at zero magnetic field and a magnetic field corresponding 
to a flux (f> — l.6(f>o with the flux quantum cf>o = hc/e. The difference marks the 
weak localization correction (from Ref. [H]). 

Fig. [2] shows the quantum transmission, numerically obtained from Eq. (J2j), 
for a "graphene billiard" [Hj, fabricated by cutting a cavity out of a two- 
dimensional graphene flake, a monoatomic layer of carbon atoms arranged in a 
honeycomb lattice. Two major quantum features are visible: (i) distinct "bal- 
listic conductance fluctuations" as a function of energy, (ii) When subject to 
an additional, perpendicular magnetic field B with magnetic flux 0, the average 
transmission (straight dashed line in the right panel) shows a small positive off- 
set compared to the average transmission for B = (solid line). This reduction 
of the average conductance at zero magnetic field reflects a weak-localization 
effect [E] . Its origin is non-classical and due to wave interference. 

Here we focus on this ballistic weak localization effect and present its semi- 
classical derivation for conductors with classically chaotic analogue. The semi- 
classical approximation enters in two steps: First, we replace G(x' ,y' ,x,y; E) 
in Eq. @ by the semiclassical Green function (in two dimensions) [3J: 

G-(r>,r;E) = £ D t{ r>, r) exp (3) 

It is given as a sum over contributions from all classical trajectories t connecting 
the two fixed points r and r' at energy E. In Eq. (J3]), 

S t (r',r;E)= [ p-dq (4) 
Jc t 

is the classical action along a path C t between r and r' and governs the accu- 
mulated phase. 

Second, we evaluate the projection integrals in Eq. ([2j for isolated trajecto- 
ries within the stationary-phase approximation. For leads with hard- wall bound- 
aries the mode wavefunctions are sinusoidal, 4> m {y) — y/2/W\ sin(mny/Wi). 
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Hence the stationary-phase condition for the y integral requires [9] 



dS_ 

dy /f , 



flT = ~Py 



"WT 



(5) 



with m = ±m. The stationary-phase solution of the y' integral yields a cor- 
responding "quantization" condition for the transverse momentum p y >. Thus 
only those paths which enter into the cavity at (x, y) with a fixed angle sin# = 
±mir/kWi and exit the cavity at (x',y') with angle sin#' = irwr/JfcWb con- 
tribute to t nm (E) with k = y/2mE/H. There is an intuitive explanation: the 
trajectories are those whose transverse wave vectors on entrance and exit match 
the wave vectors of the modes in the leads. One then obtains for the semiclas- 
sical transmission amplitudes 



'urn 

(k) 



2 A /WiWV.ftL. 

r(n,m) 



sgn(n) sgn(m) y/~At exp 



-S t (n,m; k) 



Here, the reduced actions are 

St{n,m; k) = St(k) + hky sin 6 — hky sm0 , 



(6) 



(7) 



which can be considered as Legendre transforms of the original action functional. 
The phases fit contain both the usual Morse indices and additional phases arising 
from the y, y' integrations. The prefactors are A t — \ (dy/d9') g |/ [Kk \cos9'\). 
The resulting semiclassical expression for the transmission and thereby the con- 
ductance (see Eq. ([I])) in chaotic cavities involves contributions from pairs of 
trajectories t, t'. It reads [SKUE] 



T(k) 



m— 1 n—1 



Ni N 2 



2W{W 2 



m— 1 n—1 t.t' 



with 



exp 



S t >) 



(8) 



(9) 



where the phase Ht,t' accounts for the differences of the phases jit and the sgn 
factors in Eq. ©. 

In the next two sections we show how one can calculate with purely semi- 
classical methods quantum corrections to the transmission beyond the diagonal 
approximation. The results obtained for chaotic conductors are consistent with 
RMT-predictions. 



3 Quantum transmission: 

configuration space approach 

In evaluating the off-diagonal (interference) contributions to the quantum trans- 
mission we present two approaches. The first approach is based on the analysis 
of off-diagonal pairs of trajectories and their self-intersections in configuration 
space. This approach is more illustrative, but less general than the second, 
phase space approach outlined in Sec. [U 
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We start from the semiclassical expression for \t nm (fc) | obtained by squaring 
the semiclassical transmission amplitudes, Eq. ©, see also Eqs. (I8I9|) : 



\tnm(k) | 



Tift 



2WiW 7 



^ \] A t A t > exp 



St') 



(10) 



Due to the action difference of the trajectories t and t' in the exponential in 
Eq. (I10p . this expression is a rapidly oscillating function of k or the energy in 
the semiclassical limit of large ratios (St — St') /ft- In the following we wish to 
identify those contributions to Eq. (|10p which survive an average over a classi- 
cally small but quantum mechanically large fc-window Afc. Such contributions 
must come from very similar trajectories. Afterwards we wish to evaluate their 
contributions to \t nm (k)\ 2 using basic principles of chaotic dynamics: for our 
calculation we will need hyperbolicity and ergodicity. Hyperbolicity implies 
the possible exponential separation of neighboring trajectories for long times 
with the distance growing proportional to e AT , with the Lyapunov exponent A 
and the time T of the trajectory. The second principle, ergodicity, means the 
equidistribution of long trajectories on the energy surface at the energy E of 
the trajectory. 

This section is divided into four subsections: after calculating the diagonal 
contribution we determine the simplest nondiagonal contribution in the second 
part. The behavior of the latter as a function of a magnetic field and a finite 
Ehrenfest time is then studied in the third and fourth subsection, respectively. 



3.1 Diagonal contribution 

The first, "diagonal" (D) contribution to Eq. (|10|) originates from identical tra- 
jectories t = t' ', meaning St — St'- It gives 



2W{W 2 



t 



The remaining sum over classical trajectories in Eq. (fTTj) can be calculated using 
a classical sum rule [7] that can be derived using ergodicity, see for example [16] . 
It yields 

" (T) - < 12) 

where X(-E) denotes the phase space volume of the system at energy E, and 
p(T) is the classical probability to find a particle still inside an open system 
after a time T. For long times the latter decays for a chaotic system exponen- 
tially p(T) ~ e~ T l TD , with the dwell time td — 2nh{Ph+No) ■ 1 ms exponential 
decay can be easily understood based on the equidistribution of trajectories: the 
number AN of particles leaving the system during AT is given by the overall 
number of particles N times the ratio of the phase space volume from which 
the particles leave during AT and the whole phase space volume of the system. 
The differential equation for AT, obtained in the case of infinitesimal AT, has 
obviously an exponential solution. 

By inserting Eq. (fT2"|) into Eq. (JTTJ) , we obtain 

\t nm (k)\l = w ^-. (13) 
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Figure 3: Pair of two periodic orbits in the hyperbola billiard differing from 
each other essentially in the region marked by the circle, where the left orbit 
exhibits a self-crossing while the right partner orbit shows an " avoided crossing" 
(Courtesy of M. Sieber). 

Its derivation required ergodicity valid only for long trajectories. We will assume 
that the classical dwell time is large enough, in order to have a statistically 
relevant number of long trajectories left after time t. The result in Eq. (fT3"|) 
allows for a very simple interpretation: It is just the probability of reaching one 
of the Ni + N2 channels if each of the channels can be reached equally likely 

3.2 Nondiagonal contribution 

In the following we will calculate the contributions from pairs of different trajec- 
tories, however with similar actions. For a long time it was not clear how these 
orbits could look like. There are many orbit pairs in a chaotic system having 
accidentally equal or nearly equal actions. However, in order to describe uni- 
versal features of a chaotic system after energy averaging, one has to find orbits 
that are correlated in a systematic way. These orbits were first identified and 
analyzed in 2000 in the context of spectral statistics [T7] . There, periodic orbits 
were studied to compute correlations between energy eigenvalues of quantum 
systems with classically chaotic counterpart. Based on open, lead-connecting 
trajectories, in Ref. [7] this approach was generalized to the conductance we 
study here. Still, the underlying mechanism to form pairs of classically corre- 
lated orbits is the same in the two cases. 

In Fig. [3] we show a representative example of such a correlated (periodic) 
orbit pair in the chaotic hyperbola billiard. The two partner orbits are topolog- 
ically the same up to the region marked by the circle where one orbit exhibits 
a self-crossing (left panel) while the partner orbit an 'avoided' crossing (right 
panel). Usually, such trajectory pairs are drawn schematically as shown in Fig. 

m 

One considers very long orbits with self-crossings characterized by crossing 
angles e <C it. In Ref. [T7j it was shown that there exists for each orbit a partner 
orbit starting and ending (exponentially) close to the first one. It follows the 
first orbit until the crossing, avoids this, however, traverses the loop in reversed 
direction and avoids the crossing again. 

In order to quantify the contribution of these trajectory pairs to Eq. (flQ|) 
we need two inputs: an expression for the action difference and for the density 
quantifying how often an orbit of time T exhibits a self-intersection, both quan- 
tities expressed as a function of the parameter e. The formula for the action 
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Figure 4: Schematic drawing of a pair of orbits yielding the first nondiagonal 
contribution to the transmission considered in Ref . [7J . One of the orbits crosses 
itself under an angle e, the other one possesses an "avoided crossing". Except 
for the crossing region both orbits are almost identical. 



difference AS can be derived by linearizing the dynamics of the orbit without 
crossing around the reference orbit with crossing, giving, in the limit e <C it, [17] 



AS = 



2 2 

p e 
2mA 



(14) 



At this point we can justify our assumption of small crossing angles e: In 
the limit h — > 0, we expect important contributions to Eq. (|10|) only from orbit 
pairs with small action differences, i.e. small crossing angles, as we see from Eq. 

Before deriving the number of self-crossings, P (e, T) de, in the range between 
e and e + de of an orbit of time T, we give rough arguments how this expression 
depends on e and T for trajectories in billiards. There, each orbit is composed 
of a chain of N chords connecting the reflection points. Following an orbit, the 
first two chords cannot intersect, the third chord can cross with up to one, the 
fourth chord with up to two segments, and so on. Hence, the overall number of 
self-crossings will be proportional to ^2n =3 (n — 2) cx N 2 , to leading order in N, 
i.e. proportional to T 2 . 

The crossing angle dependence of P(e,T) can be estimated for small e as 
follows: Given a trajectory chord of length L, a second chord, tilted by an angle 
e with respect to the first one, will cross it inside the billiard (with area of order 
L 2 ) only if the distance between the reflection points of the two chords at the 
boundary is smaller than Lsine. The triangle formed in the latter case includes 
a fraction sin e of the entire billiard size. From this rough estimation we expect 
P(e,T) cx T 2 sine. 

More rigorously, the quantity P (e, T) de, can be expressed for an arbitrary 
orbit 7 as [UJ 



P(e,T) de 



rT-T min (e) 



dt, 



T-t, 



,(0/2 



\JT min (e) JT min (e)/2 

xS(e-a (t s ,t s + ti))) de 



dt t \J\5(x (t.)-x(t. + fj)) 

(15) 



with the average (. . .) taken over different initial conditions (xo,po). The time 
of the closed loop of the trajectory is denoted by ti and the time before the 
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loop by t s . a(t s ,t s + U) denotes the absolute value of the angle between the 
velocities v (t s ) and v (t s + ti). | J| is the Jacobian for the transformation from 
the argument of the first delta function to ti and t s ensuring that P (e, T) de 
yields a 1 for each crossing of 7. With the absolute value of the velocity, v, it 
can be expressed as 

|J| = |v(i a ) x v(t. + i,)| =v 2 sina(t s ,t s + t l ) . (16) 

As the derivation [17] of the formula for P (e, T) for a chaotic system, starting 
from the formal expression (|15j) . is instructive to see how information can be 
extracted from the basic principles of chaotic dynamics beyond the diagonal 
approximation, we will present it here in detail. Hyperbolicity will yield a 
justification for the minimal time T m j n (e) already introduced in Eq. ()15[1 : we 
will come back to that point later and first study the effect of ergodicity. 

To proceed we interchange the phase space integral of the average with the 
time integrals, substitute (x(t s ) , p(t s )) 1 — ► (xq,Pq) in Eq. (fT5|) and obtain 

f T-T n , ID (e) 

P(e,T) = 2m dt lV 2 sme PE (e,t t ) (T - t, - T mln (e)) , (17) 

JT mm (e) 

with the averaged classical return probability density 

PB(e,tO = ^<«(xto-x(t,))<5(e-|Z(vo,v (t,))|)> ■ (18) 

This yields the probability density that a particle possessing the energy E re- 
turns after the time t\ to its starting point with the angle |Z (vo, v (t;))| = e. 
For long times this can be replaced by 1/S (E), assuming ergodicity. Then we 
obtain 

r T-T min (e) 1 
P(e,T) = 2m dt lV 2 sine— — (T - t, - T min (e)) 

= |^sine(T-2T min (e)) 2 . (19) 

Now we return to our assumption of hyperbolicity and explain the cutoff 
time T^in (e), introduced in the equations above. To this end, we consider two 
classical paths leaving their crossing with a small angle e. The initial deviation 
of their velocities is 8vi — ev. In order to form a closed loop, the deviation of the 
velocities Svf, when both paths have traversed half of the closed loop, has to be 
given by 5vt = cv with c of the order unity. Then we get for the minimal time 
Tmin (e) to form a closed loop, due to the exponential divergence of neighboring 
orbits, 

c = ee (AT mm(£ ))/2 ; (2Q) 

implying 

T min (e) = |ln(^). (21) 

An argument similar to the one used here for the closed loop can also be applied 
to the other two parts of the trajectory leaving the crossing with an angle e 
towards the opening of the conductor. Suppose t s has a length between and 
T m in(e)/2, then both parts have to be so close together that they must leave 
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both through the same lead. As we are interested in the transmission, we also 
have to exclude that cas^E A similar argument holds for the case where the last 
part of the orbit has a length between and T m j n (e)/2, in this case the orbit has 
to come very close to the opening already before the crossing and leave before 
it could have crossed. Accounting for all these restrictions, gives the integration 
limits in Eq. (|15D . 

Now we are prepared to calculate the contribution of the considered trajec- 
tory pairs to the transmission. Therefore we keep in Eq. (|ll)p one sum over 
trajectories that we will perform using the same classical sum rule as in diag- 
onal approximation, the other we can replace by a sum over all the partner 
trajectories of one trajectory, which can be calculated using P (e,T). There is, 
however, one subtlety concerning the survival probability p(T) in the sum rule: 
we argued already, that if the crossing happens near the opening, both parts of 
the orbit act in a correlated way; p(T) is changed in the case of the trajectory 
pairs considered here for a similar reason: because we know that the two parts 
of the orbit leaving the crossing on each side are very close to each other, the 
orbit can either leave the cavity during the first stretch, i.e. during the first 
time it traverses the crossing region, or cannot leave at all. This implies that 
we have to change the survival probability from p(T) to p{T — T m ; n (e) Then 
we arrive at the loop (L) contribution 



with the sum over the partner trajectories P in the first line. As the important 
contributions require very small action differences, i.e. very similar trajectories, 
and as the prefactor A t is not as sensitive as the actions to small changes of the 
trajectories, we can neglect differences between t and t' in the prefactor. In the 
second line we applied the classical sum rule with the modification explained 
before Eq. (|22|) and used P (e, T) to evaluate the sum over P. After performing 
the simple time integral in the third line, we can do the e-integration as for 
example in Ref. [18j by taking into account that the important contributions 

1 I{ we would calculate the reflection instead of the transmission, the effect of short legs, 
referred to as coherent backscattering, has to be taken into account. 

2 This effect together with the requirement of a finite length of the orbit parts leaving 
towards the opening was originally not taken into account in Ref. [7]. In this calculation, 
the contributions from these two effects cancel each other, they will be only important when 
considering more complicated diagrams as in the next section. 





(22) 
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come from very small e, yielding 



/,m2 8nhmv 2 T% C , , i \ ~ ■ ( P 2 £ 2 

\tnm(k)\ L = — — / de(e/c)^D sine cos 



£(£) 2 J w ' \2m\h 
8irhmv 2 Tn f , m\h ( 1\( T7 d) / 2m\hz\' 



J2. I dz 

£(£) 2 J p 2 \c 



8nhmv 2 Tn f , rah I 1\ \ Xt d) ( 2mXHz\ Xt o smz , , 

- -^f/^U (—J — ' (23) 

In the first line we already rewrote e _Tmin ^^ ri3 as (e/c) Xt ° , and in the second 
line we approximated sine w e and substituted 2 = p 2 e 2 / (2m\h). Then we 
perform a partial integration with respect to z neglecting rapidly oscillating 
terms that are cancelled by the fc-average, introduced after Eq. (fTTJ|) . Eventually, 
we perform the z-integral by pushing the upper limit to infinity, i.e. h — > 
and taking into account our assumption of large dwell times, i.e. Xtd — > 00. 

Additionally we assume (2mAfi/p 2 ) ~ 1; we will return to the last point 
soon. 

Finally, we arrive at the leading nondiagonal contribution to the quantum 
transmission [7], 

\t nm (k)\ 2 L = -- -2. (24) 

(N1 + N2) 

3.3 Magnetic field dependence of the nondiagonal contri- 
bution 

Up to now we assumed time reversal symmetry. If this symmetry is destroyed, 
for example by applying a strong magnetic field, the latter contribution will 
vanish, because the closed loop has to be traversed in different directions by the 
trajectory and its partner. Here we study the transition region between zero and 
finite magnetic field. In particular, we assume a homogeneous magnetic field 
B z perpendicular to the sample that is assumed weak enough not to change the 
classical trajectories, but only the actions in the exponents. Since the closed 
loop is traversed in different directions by the two trajectories, we obtain an 
additional phase difference (AttAB z / '(f>o) between the two trajectories with the 
enclosed area A of the loop and the flux quantum cf>o = (hc/e). We further need 
the distribution of enclosed areas for a trajectory with a closed loop of time T 
in chaotic systems, given by 

F ( AT > = 7m exp (-^) (25) 

with a system specific parameter (3. A derivation of this formula can be found 
for example in Ref . |19j . Including the phase difference and the area distribution 
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in a modified P(e,T) yields 



Pb (e, T) = — — sin e / dfc, (T — tj — T min (e)) 

/OO 
dAP(^,t;-T min (e))cos- 
-OC 



2m V 2 /^-We) 
■ sin e 



£(J5) 



dti (T- t, - T min (e)) e -(* ! -^(e))/t B(26 ) 

We) 



with ts = S7r fp B 2 ■ In the first line we used that paths leaving the crossing to 
form a closed loop, enclose a negligible flux, as long as they are correlated; for 
a more detailed analysis see Appendix D of Ref. [30]. Performing the T- and 
e-integrals similar to the case without magnetic field, yields [7] 

\t nm (k,B z )\ 2 L = - 1 1 . (27) 

(iVi + iVa) 1 + T D/tB 

We obtain an inverted Lorentzian with minimum at zero magnetic field, implying 
that the transmission through our sample increases with increasing magnetic 
field. This weak localization phenomenon, a precursor of strong localization, is 
visible as the reduction of the average quantum transmission in Fig. [2J 

3.4 Ehrenfest time dependence of the nondiagonal contri- 
bution 

This semiclassical approach can also be applied to calculate the Ehrenfest time 
dependence of the transmission. The Ehrenfest time te = In (E/(\K)), 

more generally a time proportional to In ft |21j . is the time a wave-packet needs 
to reach a size such that it can no longer be described by a single classical par- 
ticle. The Ehrenfest time thus separates the evolution of wave packets following 
essentially the classical dynamics (e.g. up to few bounces of the wave packet in 
Fig. [T]) from longer time scales dominated by wave interference (last panel in 
Fig. QJ. Based on field theoretical methods, Aleiner and Larkin showed, that a 
minimal time is required for quantum effects in the transmission to appear, the 
Ehrenfest time 22J. We will now apply our semiclassical methods to determine 
the Ehrenfest-time dependence of the transmission following the pioneering work 
[2"5] that has later been extended to the reflection O [551 HO] , including a dis- 
tinction between different Ehrenfest times. To this end, we directly start from 

Eq. (|22p ; however, we have to be more careful when evaluating the e- integral: 

, i 

in our former calculation we assumed e" TE / TD = (\h/E)^ D w 1, requiring 
te <§C Try. Lifting however this strong restriction for td, but still keeping it 
large enough to fulfill our assumption of chaotic dynamics, we obtain 

\t nm {k,T E )\l = - — 1 AT , 2 e- TE/TD , (28) 

i.e. an exponential suppression of the nondiagonal contribution due to the 
Ehrenfest time. This dependence has been confirmed in numerical simulations 
for quantum maps. 
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After this introduction into semiclassical methods for the evaluation of non- 
diagonal contributions in configuration space, we will now turn to the generaliza- 
tion to phase space which also allows for an elegant way to compute higher-order 
corrections to the leading weak localization contribution presented above. 

4 Quantum transmission: 
phase space approach 

The above configuration space treatment, based on self-crossings, is restriced 
to systems with two degrees of freedom. More generally, for higher-dimensional 
dynamical systems one cannot assume to find a one-to-one correspondence be- 
tween partner orbits and crossings of an orbit 26J. In order to overcome these 
difficulties, a phase space approach was developed for calculating the spectral 
form factor in spectral statistics in [Ml H7J involving periodic orbits. The next 
challenge was the generalization of this theory to trajectory pairs differing from 
each other at several places, solved again first for the spectral form factor [5H] 
and generalized to the transport situation considered here in Ref. which 
serves as the basis of the following discussion. 

In this section we first explain the phase space approach and use it after- 
wards in the way developed in Refs. [28j ES] for the calculation of the quantum 
transmission involving also higher-order semiclassical diagrams. 

4.1 Phase space approach 

Compared to the last section, we first of all have to replace the role of the 
reference orbit. Whereas we used there the crossing orbit as reference and 
calculated then the action difference and the crossing angle distribution in terms 
of the crossing angle e, we will consider here the orbit without the crossing that 
is close to itself in the region, where the partner orbit crossed itself in the last 
section. Inside the region, where the orbit is close to itself; we will refer to it 
as encounter region and to the parts of the orbit inside as encounter stretches, 
that are connected by so called links; we place a so called Poincare surface 
of section V with its origin at x = (x, p). The section consists of all points 
x+Sy = (x + <5x, p + <5p) in the same energy shell as the reference point with the 
<5x perpendicular to the momentum p of the trajectory. For the two-dimensional 
systems considered here V is a two-dimensional surface, where every vector by 
can be expressed in terms of the stable direction e s (x) and the unstable one 
e u (2) pS] 

Sy — se s (x) + ue u (x) . (29) 

The expressions stable and unstable refer to the following: we consider two 
orbits, one starting at x and the other one at x+dyo, then the difference between 
the stable coordinates will decrease exponentially for positive times and increase 
for negative time exponentially in the limit of long time T, unstable coordinates 
behave just the other way round. The functional form of the exponentials can 
be determined as e AT and e~ AT . Now we can come back to the trajectory 
with the encounter region, where we have put our Poincare surface of section. 
The trajectory considered in the last section will pierce through the Poincare 
section twice: one of this points we will consider as the origin of the section 



13 



the other piercing will take place at the distance (s, u). The coordinates of the 
piercing points of the partner trajectory are determined in the following way: 
the unstable coordinate of the partner trajectory has to be the same as the one 
of that part of the first trajectory that the second will follow for positive times. 
The stable coordinate is determined by the same requirement for negative times. 

After this introduction into the determination of the (s, it)-coordinates, we 
are now ready to treat trajectory pairs that differ in encounters of arbitrary 
complexity. Following Ref. [29j and using the notation introduced there we will 
allow here that the two trajectories differ in arbitrarily many encounters in- 
volving an arbitrary number of stretches. In order to organize this encounter 
structure we introduce a vector v = (i>2, ^3, ■•■) with the component vi determin- 
ing the number of encounters with I stretches involved. The overall number of 
encounters during an orbit will be denoted by V — Y^b=2 Vl > tne overau number 
of encounter stretches by L = X);=2 ^ v i- ^ n an encounter of I stretches, we will 
get I — 1 (s, u)-coordinates. 

Now we can proceed by replacing the former expressions for the minimal 
loop time, the action difference and the crossing angle distribution depending 
on the former small parameter e, by the corresponding expressions depending 
on the new small parameters (s,u). 

We start with the minimal loop time, to that one also refers as the duration 
of the encounter. Shifting the Poincare surface of section through our encounter, 
the stable components will asymptotically decrease, the unstable ones will in- 
crease for increasing time. We then claim that both components have to be 
smaller than a classical constant c, its exact value will again, as in the last 
section, be unimportant for our final results. We finally obtain the encounter 
duration i cnc as the sum of the times t u , that the trajectory needs from V till 
the point where the first unstable component reaches c, and the time t s , that 
the trajectory needs from V till the point where the last stable component falls 
below c. Thus we get 

1 c 2 

tene = ts tu ~~ T" In f] Tj f"j TT- (30) 

A maxj {Isiljmaxj {\Uj\\ 

Now we treat the action difference between the two trajectories. Expressing 
the actions of the paired trajectories as the line integral of the momentum along 
the trajectory, we can expand one action around the other and express the 
result in terms of the (s, u)-coordinates. This yields for an ^-encountei|3 

i-i 

AS = J2s jUj , (31) 

j=i 

The action difference of a trajectory pair is then obtained by adding the differ- 
ences resulting from all encounters. 

Finally we come to the crossing angle distribution that will be replaced here 
by a weight function for the stable and unstable coordinates for a trajectory of 
time T. We first notice that the uniformity of the trajectory distribution implies 
in terms of our coordinates in V that a trajectory pierces through the section 



3 Strictly speaking 28] the (s, u) coordinates used here and in the following calculation are 
for encounters involving more than two stretches not the same like the ones described before, 
but related to them via a linear and volume preserving transformation. 
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with the coordinates (u, u + du) and (s, s + ds) within the time interval (t, t + dt) 
with the probability l/Y>(E)dsdudt. In general, we obtain for an Z-encounter 
(1/£(.E)) 1 ds l ~ 1 du l ~ 1 dt l ~ 1 . Integrating the product of the latter quantities 
for all encounters over all possible durations of the L — V intra-encounter links, 
in a way that their durations are positive yields our weight function for a fixed 
position of V. To take into account all possible positions of V, we also integrate 
over all possible positions, where it can be placed and divide by t cnc to avoid 
overcounting of equivalent positions. Taking all link times positive, we obtain 
for the weight function for an orbit of time T 

1 

(£(i^-~riLi*?„c 



w T (s, u) = — ,^^ L _ V — / dh . . . dt L Q \ T - 

.1 ^er 

L 




L!(E(J5)) £ - v IlLi*S« 



(32) 



One additional problem arises, when treating trajectory pairs differing not 
only in one 2-encounter like in the last section: one can construct for one v 
different trajectory pairs, varying for example in the relative orientation, in 
which the encounter stretches are traversed. We will count this number by a 
function N (v) and describe briefly later, how it can be calculated. 



4.2 Calculation of the full transmission 

After the introduction into the phase space approach we are now ready to cal- 
culate the transmission. Taking the weight function, the action difference and 
the number of structures, we can transform the non-diagonal (ND) part of Eq. 
(fTOf into 

\t nm mi D = ^EE^w 

t V 



... J" d L - v ud L - v s exp Q A s) w t (s, u)^ (33) 



Ak 

with the average over a small fc-window denoted by (. . .) Afe . Inserting the formu- 
las for the action difference, the weight function and using the classical sum rule 
with the modification of the survival probability, discussed in the last section, 
we can then transform the integral with respect to the length of the trajectory 
into one over the last link and obtain 



\tnm (fc) | 



\ND 




d L - v ud L - v s ex P (-% + 1 E-=i X ■ 



(E(£)) L - V ^ t c Q „c 



I 1 \ ST ri\V 



iVi + N 2 ^ \ N X + J ^ 

n— 1 v ' v 



L-V=n 

(-l) V N(v) (34) 
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with the L + l link times The (s, u)-integrals are calculated using the rule 
[25] . that after expanding the exponential e tcnc / Tr> into a Taylor series, only the 
t cnc -independent term contributes and yields in leading order in h 

mIJj sduexp &)) A r^ (35) 

with the so called Heisenberg time Tu — S(-B) / (2-kK). For the sum with respect 
to v, one can derive recursion relations, yielding |28j 

L-V=n , o\™ 

(-l) V N(v)=(l--\ (36) 

V 

with (3 = 1 and (3 = 2 for the case with and without time reversal symmetry, 
respectively. Relations of this kind are derived by describing our trajectories 
by permutations expressing the connections inside and between encounters and 
considering the effect of shrinking one link in an arbitrarily complicated struc- 
ture to zero. 

We then obtain for T(Ep), given in Eq. (JTJ) , in the case with time reversal 
symmetry [33] 

T (E )^ ~ NlN2 + NlN2 V ( 1 \ n = NlN2 (37) 
K F} ~ Ni+N 2 n 1 +N 2 ^\N 1 +N 2 ) Ni+N 2 + 1 1 ' 

and in the case without time reversal symmetry 

which agrees with the diagonal contribution, already obtained in the previous 
section. Both results are in agreement with RMT predictions [3T] , 



5 Semiclassical research paths: 
present and future 

In the above sections we outlined the recently developed semiclassical techniques 
for treating off-diagonal contributions for one paradigmatic example of coherent 
transport, weak localization. Nonlinear (hyperbolic) dynamics as a prerequisite 
for the formation of correlated orbit pairs, together with classical ergodicity are 
at the core of the universality in the interference contribution to the ballistic 
conductance. Recently, this semiclassical approach to quantum transport has 
been extended along various paths, thereby gaining a more and more closed 
semiclassical framework of mesoscopic quantum effects: 

In the presence of spin-orbit interaction, spin relaxation in confined bal- 
listic systems and weak antilocalization, the enhancement of the conductance 
for systems obeying time-reversal symmetry (i.e. at zero magnetic field), has 
been semiclassically predicted [32] and generalized to other types of spin-orbit 
interaction 33J and higher-order contributions |34j . 

Ballistic conductance fluctuations, the analogue of the universal conductance 
fluctuations in the diffusive case, require the computation of semiclassical dia- 
grams involving four trajectories. Again, RMT predictions could be confirmed 
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[25] . Accordingly, RMT predictions for shot noise in Landauer transport agree 
with recent semiclassical results [35J. This could also be shown for higher mo- 
ments of the conductance in [36] . Even strong localization effects (in certain 
systems, i.e. chains of chaotic cavities) could be derived by making use of semi- 
classical loop-contributions [37] , 

Finally, beyond RMT, Ehrenfest time effects on transport properties have 
recently been in the focus of the field. As a result, ballistic conductance fluctu- 
ations turn out to be, to leading order, TE-independent [24] , contrary to weak 
localization which is suppressed with e~ TE ^ TD [23 , as we have shown. Recently 
these results were extended to the ac-conductance in [38] , 

Beside universal spectral statistics of closed systems and transport through 
open systems, the semiclassical techniques have been recently generalized to 
decay and photo fragmentation of complex systems [39j . 

To summarize, by now there exists a newly developed semiclassical machin- 
ery to compute systematically quantum coherence effects for quantities being 
composed of products of Green functions. This opens up various possible future 
research directions: With respect to system classes, mesoscopic transport theory 
has been predominantly focused on electron transport, leaving the conductance 
of ballistic hole systems aside, which plays experimentally an important role and 
for which a semiclassical theory remains to be developed. The situation is simi- 
lar for a novel class of ballistic conductors: graphene based nanostructures, see 
Fig. [2 More generally, a proper treatment of interaction effects in mesoscopic 
physics would probably represent the major challenge to future semiclassical 
theory. 

To conclude, at mesoscopic scales where linear quantum evolution meets 
non-linear classical dynamics, further interesting phenomena can be expected 
to emerge, in view of increasingly controllable future experiments. 
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